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Abstract The specific heat capacity of a two-dimensional electron gas is de- 
rived for two types of the density of states, namely, the Dirac delta function 
spectrum and that based on a Gaussian function. For the first time, a closed 
i-Q ■ form expression of the specific heat for each case is obtained at half-filling. 

When the chemical potential is temperature-independent, the temperature 
is calculated at which the specific heat is a maximum. Here the effects of the 
broadening of the Landau levels are distinguished from those of the different 
filling factors. In general, the results derived herein hold for any thermody- 
Cd ' namic system having similar resonant states. 
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_ ., ' 1 Introduction 
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fvi ■ Electrons in a two-dimensional layer continue to be of interest in physics since 

they exhibit nonclassical behavior and are readily realizable in semiconductor 
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heterostructures [T]. Nowadays research investigations on two-dimensional 
electron gas (2DEG) systems are aimed at understanding the plateaus in 
the quantum Hall conductance [2] and the de Haas - van Alphen effect or 
the magnetization oscillations J3,3»S] ■ The smooth drop of the plateaus and 
the extrema of the oscillations reveal broadened energy levels, that is, a 
finite density of states (DOS) [5] between eigenvalues E n . These E n , known 
as Landau levels, are the quantized energy spectrum (E) obtained when a 
strong magnetic field (B) is applied perpendicular to the system's plane [7] or 
in a tilted direction relative to it [8]. The striking behavior of E n broadening, 
from which other exotic features arise, has been attributed to disorder due to 
the presence of impurities, defects and other inhomogeneities in the system 

mm]. 

Ideally, the DOS of a noninteracting 2DEG in a perpendicular B is given 
as a series of delta functions [7], that is, 

n 

Here e is the electron charge, h is Planck's constant and E n = (n + l/2)hcu c is 
the nth Landau energy level. The latter depends on B through the cyclotron 
frequency u c = eB/m* for a given effective mass m*. Such a DOS structure 
can be used to model actual materials with a very narrow distribution of 
energy carriers [3J. 

Usually, the actual shape of the density of states of a 2DEG is determined 
by making theoretical fits to the heat capacity data from experimental mea- 
surements. This procedure made on a 2DEG in GaAs-GaAlAs multilayers 
[T2] . for instance, yields a DOS with Gaussian peaks on a flat background. 
Many experimental and theoretical studies on the thermodynamic and mag- 
netic properties of 2DEG systems, therefore, use a form of DOS based on a 
Gaussian function [Tl lllllTSirB] such as 
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where the broadening is taken into account by the parameter r. 

In one study using this DOS, the authors of Ref. [T4] predicted that the 
specific heat and other magneto-thermal effects at high B become universal 
functions of ksT/ r, where fc^ is the Boltzmann's constant and T is the tem- 
perature. However, their numerical calculations at constant electron concen- 
tration were unable to provide the exact dependence of these thermodynamic 
properties on UbT '/ T '. 

Another study of an impurity-free quantum Hall system in a periodic 
rectangular geometry found that the specific heat at constant volume behaves 
as Cy ex exp[— A/T]/T 2 , where A is the energy gap between the ground state 
and the first excited state [TB] . The numerical calculations in Ref. [TB] derived 
Cy for different values of the filling factor v. This nominal number of filled 
Landau levels is given as v = hN/Be, where N is the electron concentration. 
For various filling factors, their results yield a Cy with a sharp peak at low 
T whose height depends on A. However, their Cy calculations show that 



finite size effects become dominant at low T. Without an exact expression 
for Cy, it would also be difficult to distinguish the exp[— A/T]/T 2 behavior 
from the dominant lattice contribution ~ T 3 at very low T. In general, since 
Cy depends strongly on the nature of a given system such as the DOS, the 
exact behavior of Cy as a consequence of its electronic structure is not easily 
determined. 

In this work, the authors will show that a closed form expression of Cy 
is obtained for a 2DEG system when the chemical potential \x does not vary 
with T. The latter condition is satisfied at a finite non-zero T when the 
last occupied Landau level is half-filled, that is, the last occupied energy 
level is filled to half its degeneracy. Particular interest has been devoted to 
v = 1/2 and other half- filled states since no plateau of the Hall resistance 
is observed unlike other fractional v with an odd denominator 17 . At half- 
filling, the chemical potential is right at the center of the Landau level where 
extended states are located [2J. There is an absence of a gap at this v and 
some experiments have observed features of a Fermi-liquid state [TSinU] . 

Here a detailed derivation will be presented of Cy at v = n'/2 (where n' 
is an odd integer) for a spinless 2DEG having a delta-shaped or a Gaussian 
DOS in the presence of a perpendicular applied magnetic field. The exact 
relations for the specific heat as obtained for the first time in this work serve 
as a basis that characterizes the underlying density of states of an electron 
gas. In addition, the condition at which the heat capacity is a maximum is 
established for both cases including the actual nature of its dependence on 
r. This will be particularly useful for the design of devices that are intended 
to resist changes in temperature. Lastly, remarks will also be given regarding 
Cy at other filling factors and hence, the general case will be formulated. 



2 The Temperature Dependence of the Chemical Potential 

The chemical potential fi, which measures the amount of energy required to 
add a particle in a given system, is generally a function of T and B. But at 
half-filling the /i is independent of T since the 2DEG is metallic and there 
is no energy gap between the conduction and the valence bands. This is 
demonstrated in Fig. [T] The chemical potential shown here is obtained for a 
fixed electron concentration, 

N = J f(E)D(E) dE , (3) 

using the root- finding method as outlined in Ref. [15]. Here f(E) is the 
Fermi-Dirac distribution function and it is given as 

f{E) = exp[(£ - fi)/k B T} + 1 • (4) 

The value of N is set equal to 3.6 x 10 11 cm -2 , which has the same or- 
der of magnitude as in experiments of 2DEG systems in semiconductors. In 
addition, the 2DEG system used has a Gaussian DOS as given in Eq. ([2]) 
with r — 0.2 meV. No spin-splitting is considered here. Hence, in Fig. [T] 
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Fig. 1 The temperature-dependent chemical potential at various filling factors for 
a 2DEG system with a Gaussian DOS. Here u* — 5 corresponds to half-filling for 
spinless electrons. This illustrates that fi does not change with T in this case. As 
examples, the inset provides the derivative of /j, with respect to T for v* — 4.5 and 
v* =6. 



the even integer v* corresponds to full filling and the odd integer v* to half- 
filling, that is, v* /2 — > v. As the T approaches zero, the /i of the fully-filled 
states (v* = 4, 6) varies within huj c /2 from the temperature-independent 
\x (y* — 5). This is likewise observed for the non-half partially filled states 
(u* = 4.5, 5.5), albeit at a gradual slope. This signifies that the [i behavior 
and, consequently, the specific heat nature of 2DEG systems at the non-half- 
filled states can be considered as deviations from their behavior at half-filling. 
In fact, the /j, saturates to its half-filling value as T increases. Treating the 
fj, as T-independent is a valid approximation even for non-half-fillings albeit 
limited only at high T. 



3 The Heat Capacity at Half-Filling 

The specific heat of a substance is a measure of the amount of energy needed 
to raise its temperature by a degree. For an electron gas at constant volume, 
Cy can be obtained from its internal energy U as follows 



C v = 



dT ~ dT 
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When the ^ does not vary with T, only the f(E) remains as the In- 
dependent function. The temperature derivative in Eq. ([5]) then acts only on 



f(E) and Eq. ([5]) can be rewritten as 
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3.1 For a Delta-shaped DOS 

As shown in the numerical calculations of Ref. [15] at constant N, despite 
the presence of the broadening of Landau levels, a characteristic low tem- 
perature is reached at which the 2DEG approaches its ideal electron gas 
behavior. Therefore, the continued use of a delta-shaped DOS is justified in 
heat capacity measurements at strong B and below a characteristic low T. 

Substituting D(E) of Eq. fTJ) in Eq. fl£J), we obtain the specific heat at 
half-filling for a delta-shaped DOS, that is, 

B h ^cosh 2 [(E n -^)/2k B T] 

Note that the number of occupied Landau levels and the chemical potential 
are determined if the electron concentration is known. For this case, the 
contribution of the last occupied Landau level, n — n max , to the CV is zero 
since E n = /i. 

The specific heat expressed in Eq. ([7} exhibits a single peak at the low 
temperature region. The temperature at which CV is a maximum can be 
obtained by taking dCy/dT = and performing a second-derivative test. 
The non-trivial condition at which this is satisfied is when 



E< 



^ tanhK) ' [K) _ coth((Jw)]=0> (8) 

cosh (uj n ) 

where w n = (E n —/j,)/2kBT and E n =/= /i. The maximum specific heat becomes 

Cv ^^ kB — g sinh 2 [(£„- M )/2fc B Tf (9) 

The value of T when Cy is a maximum, TUg^, can be determined at 
n = n max — 1 because this is the most dominant term. Also, i?„ mai _i ^ fi 
and Eq. ((5]) remains valid. At this level, solving for ui n ~ coth(w„) = yields 
u n = 1.19968. Thus, we find that T peak ~ 0.416778(K mo3: -i - n)/k B - This 
result also provides a measure of the Landau level spacing. 

The derivations herein are valid for any B and for any fi(B) as long as 
the condition of half-filling is met or whenever the chemical potential does 
not change with temperature. 



3.2 For a Gaussian-shaped DOS 



There are other forms of the density of states used in the literature [3l|4j 
I12j such as the Lorentzian, Gaussian with a constant background and the 
semielliptical distribution function. Here we consider only the Gaussian DOS 
as given in Eq. ^ to demonstrate the effect of a finite broadening r as 
compared to the ideal impurity-free case. We keep in mind that in the limit 
as the r — > 0, the Gaussian DOS approaches the delta-shaped DOS as given 
in Eq. (JXJ) since 

5{x) = lim r ^o (l/2 7 rr 2 ) 1 / 2 exp(-x 2 /2r 2 ). 
Substituting Eq. © in Eq. © yields 
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(10) 



For clarity and brevity, dimensionless energy parameters lo — (E—fi)/2k B T 
and uj n — (E n — ^i)/2k B T are introduced such that the specific heat given in 
Eq. (fTU|) reduces to 
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Here 7 = (l/2) 1 / 2 7^/fcsT is the dimensionless broadening parameter. 

The integral in the specific heat expression, Eq. (jlip . is evaluated in two 
sequences of integration by parts which involve the integrals from Rcf. 20 : 
that is, (i) following Integral No. 8 of Section 2.477 [20] 
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and that of (ii) Integral No. 7 of Section 2.479 [20] 
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where \u\ < tt/2 andp > — 1. The constant B 2 k is the 2fc-th Bernoulli number. 
The infinite integration of products of powers and hyperbolic functions are 
best treated using the Bernoulli series as a generating function in regions 
where the hyperbolic functions converge. 

In terms of these integrals, the specific heat for a Gaussian DOS is given 
as 
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The remaining integrals in Eq. (J14I) are of the form 
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where for the case when j is odd 
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and for the case when j is even 
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In this equation, erf[z] = (2/V7T) J exp[— z 2 } dz is the error function. The 
series and integrals of Eq. (I14p converge. Therefore, one can determine the 
exact value of CV at a given T and B for a particular 7. 

Note that, at half- filling, E n = /i at the last occupied Landau level, 
n = n max . In these cases, w n = 0. Hence, the integral of Eq. (IT51) van- 
ishes. Therefore, the specific heat at the last occupied Landau level given by 
Eq. (fT4|) evaluated from cj = — tt/2 tow = tt/2 at cj„ = yields 
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where the coefficients C\ and C 2 are constants. From Eq. (|T4|) . these constant 
coefficients are derived to be as foilows 
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Their numerical values are C\ — 5.07023 and C 2 — 3.6412. 

Recall that 7 is inversely proportional to T. Setting to zero the derivative 
with respect to T of Eq. ([T5]) . a quartic equation is obtained. This yields one 
real root with T > 0. This is the temperature at which the specific heat has 
a maximum at E n — fi, namely, 
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where C 3 = vr 4 tanh[7r/2] = 89.339. This yields that T peak = 0.501815 r/k B . 
It is important to focus on the last occupied Landau level at n = n maa; 
at high B and low T because this is where the significant contribution to Cy 
comes from [13] ■ Since we are considering half- filling, E nmax = fx, the specific 
heat as derived from Eq. (ITU1) becomes Eq. ([T5|. When this is rewritten in 
terms of k B T / r, we find that 
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where the coefficients C\ and C2 are constants. They are given in Eq. (|19p and 
Eq. ([2"U|) . respectively. This result is consistent with the case for bulk metals 
that the Cy ~ aT + 6T 3 where a and 6 are arbitrary coefficients. The linear 
in T behavior of Cy comes from the electronic contribution, while the T 3 
behavior is associated with lattice vibrations. This broadening of the Landau 
levels as represented by a Gaussian DOS that yields the T 3 contribution to 
the specific heat may be in part of the same nature as phonons. In addition to 
the similarity with metallic bulk properties, we obtain a Gaussian dependence 
on k B T/r for the specific heat as a natural consequence of the DOS, that is, 
CV - exp[-(fc B T/r) 2 ]. This is unlike in Ref. [Tg] where C v ~ exp[-A/T}. 
Thus, in this work, the drop in Cy after the T pea k is caused by the Gaussian 
tail. 

Equation (|22l) serves as a reference for the Cy measurement as T is varied. 
When r is comparable to or greater than k B T, Cy decreases as r increases 
for a given T. This is due to the lesser energy requirement to move from one 



energy level to the next since there is a larger overlap between neighboring 
Landau levels. When T « k B T, the exponential behavior of Eq. (l2"2l dom- 
inates. In this case, Cy increases as T increases. We have assumed here that 
r is a constant for a given B and T. Also, T is assumed to be the same for all 
Landau levels. However, the effect of T will be balanced by an increasing B 
that increases the separation between Landau levels while raising the degen- 
eracy of the states. Even when a -B-dependent critical value of r is reached 
wherein the blend of states disappear [5] , the expression for Cy t „ max will still 
hold. 

Furthermore, the functional form of Cv in Eq. (|2"2")l yields a single peak. 
As obtained from Eq. (|2ip . this occurs at T pea k = 0.501815 T/k B . This linear 
proportionality with T Xiea ^ L and Eq. (|22p can be used to directly measure the 
broadening parameter of the DOS for a given B whenever experimental data 
for the CV in a 2DEG system are available. 



4 A Note on the General Case 



When the chemical potential is temperature-dependent, the derivative of the 
Fermi function is of the general form 



(23) 



This will yield two terms in the specific heat as given in Refs. [T3"UT4] . But 
aside from f(E), the specific heat also has an explicit dependence on d/j,/dT. 
Evaluating the temperature derivative in Eq. ([5]) results into 
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Using in this expression the Eq. ((23)) , we obtain the detailed formula for Cy 
when [i — (a(T) as follows 



Cy — k 
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This is the general equation for the specific heat and is valid regardless of 
the filling factor used. Hence, we find that the correction to Eq. ([5]) for 
Cy will depend on how large d^/dT is. In the case of a Gaussian-shaped 
DOS, we have numerically shown that the contribution from d\ijdT for a 
fixed concentration is small. In the example given in the inset of Fig. [I] this 
contribution is of the order of 0.01 — 0.1 meV/K. 

In the case of the delta-shaped density of states, substituting the DOS of 
Eq. (P, in Eq. © yields 
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where u> n = (E n — /x)/2fcgT. We get two corrections terms to Eq. Q. Their 
values depend again on the magnitude of d/i/dT. 

Though the focus of this work is on the 2DEG, the theoretical model as 
presented is also applicable for any system under the grand canonical ensem- 
ble formulation with either a delta-shape or a Gaussian-broadened density 
of states. 



5 Conclusions 

The cardinal behavior of the specific heat of two-dimensional electron gas 
systems in an applied perpendicular magnetic field is derived for the cases 
wherein the density of states either has a delta-shape or a Gaussian form. 
This work establishes for both cases of the DOS the leading contributions 
to the specific heat and the temperature condition at which the Cy is a 
maximum. The exact expressions obtained for Cy hold when the chemical 
potential is independent of the temperature or the 2DEG is in a half-filling 
state. At low T and when the fi depends on T, the formulation derived 
herein remains valid with an additional term in Cy that is proportional to 
the d[ijdT . Furthermore, the results provide the linear proportionality of the 
broadening of the Landau levels to the temperature at which the Cy attains 
its highest value. Experimental realizations of the results for the broaden- 
ing parameter can be made, for example, if it is induced by the presence 
of impurities, by varying impurity levels of a half-filled 2DEG system and 
monitoring the corresponding changes of the temperature at which the Cy is 
a maximum. When effects of various spin configurations are then considered, 
the derivations in this work can be used to distinguish them from the effects 
of the temperature dependence of the /j, and the broadening of the Landau 
levels. 
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